Theoretical Study of Carbon Clusters in Silicon Carbide Nanowires 
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Using first-principles methods we performed a theoretical study of carbon clusters in silicon car- 
bide (SiC) nanowires. We examined small clusters with carbon interstitials and antisites in hydrogen- 
passivated SiC nanowires growth along the [100] and [111] directions. The formation energies of 
these clusters were calculated as a function of the carbon concentration. We verified that the en- 
ergetic stability of the carbon defects in SiC nanowires depends strongly on the composition of the 
nanowire surface: the energetically most favorable configuration in carbon-coated [100] SiC nanowire 
is not expected to occur in silicon-coated [100] SiC nanowire. The binding energies of some aggre- 
gates were also obtained, and they indicate that the formation of carbon clusters in SiC nanowires 
is energetically favored. 

PACS numbers: 71.15.Nc, 73.22.-f, 61.72.J- 
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I. INTRODUCTION 

Silicon carbide (SiC) is a wide-band-gap semicon- 
ductor with excellent physical, electronic and mechani- 
cal properties^ such as high thermal conductivity, high 
breakdown field, low density, high saturation velocity, 
high mechanical strength, and stability at high temper- 
ature. These exceptional features make SiC a promising 
candidate to replace silicon in electronic devices operat- 
ing in high-power, high-frequency, and high-temperature 
regimes^ 

In the last years, SiC nanostructures (like 
nanospheresjS nanospringSft nanowires 5 - and nanotubes^) 
have been successfully synthesized, and several theoret- 
ical and experimental works 3 - - — have been performed 
to investigate their structural and electronic properties. 
The unique features of SiC combined with quantum-size 
effects make the SiC nanostructures interesting ma- 
terials for nanotechnology applications. For instance, 
SiC nanowires and nanotubes have been considered as 
candidates for hydrogen storage nanodevices^ and for 
building blocks in molecular electronic applications. 12 
In particular, silicon carbide nanowires (SiC NWs) have 
excellent field emission properties^ high mechanical 
stability and high electrical conductance* 7 - and could be 
used as nanoscale field emitters or nanocontacts in harsh 
environments. 

Some optical and electronic properties of semiconduc- 
tors may be modified by the presence of defects. The 
most common defects in SiC are vacancies, intersti- 
tials, antisites and clusters. These defects are mainly 
formed during the growth process and ion implantation 
of dopants. Vacancies and interstitials of C and Si in 3C-, 
4H- and 6H-SiC bulks have been thoroughly investigated 
in theoretical and experimental worksJ^— These inves- 
tigations have showed that the C and Si vacancies are 
electron and hole traps,— - — and that C and Si intersti- 
tials have higher mobility than vacancies, although the 
mobility of point defects in SiC is reduced as compared to 



another semiconductors (like silicon) The high mobil- 
ity of carbon interstitials favors the formation of carbon- 
interstitial clusters. Using ab initio methods, Gali et al. 20 
systematically investigated small clusters of carbon in- 
terstitials and antisites in 3C- and 4H-SiC bulks, and 
verified that the formation of carbon aggregates is ener- 
getically favored. 

In spite of some theoretical studies on carbon aggre- 
gates in SiC bulk^ - — the investigation of carbon clus- 
ters in SiC NWs is very scarce. Motivated by the lack of 
studies on C aggregates in SiC NWs and by their ener- 
getically favorable formation in SiC bulk, in this work we 
performed an ab initio study of small carbon clusters in 
SiC NWs. We considered hydrogen-passivated SiC NWs 
grown along the [111] and [100] directions and examined 
clusters with interstitial and antisites carbon atoms. The 
formation energies of these clusters were determined as a 
function of the C concentration. We calculated the bind- 
ing energies of some aggregates, and our results indicate 
that the formation of carbon clusters in SiC nanowires is 
energetically favored. Besides [111] and [100] SiC NWs, 
the carbon clusters were also investigated in 3C-SiC bulk, 
in order to compare the effect of C defects in SiC bulk 
and NW. 



II. METHODOLOGY 

In this work we present ab initio calculations 
based on density functional theory 2 * 3 (DFT) car- 
ried out by using the SIESTA code^ We used lo- 
cal spin-density approximation 25-27 (LSDA) for the 
exchange-correlation functional and norm-conserving 
fully-separable pseudopotentials 2 ^ to treat the electron- 
ion interactions. The Kohn-Sham orbitals were expanded 
using a linear combination of numerical pseudoatomic 
orbitals 2 ^. and a double-zeta basis set with polarization 
functions 3 ^ (DZP) was employed to describe the valence 
electrons. 
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The 3C-SiC bulk, and the [111] and [100] SiC NWs, 
were modeled within the supercell approach, with 128, 
232 and 279 atoms, respectively. The SiC NWs were 
constructed from the 3C-SiC structure and the dangling 
bonds of their surfaces were saturated with hydrogen 
atoms. Due to the periodic boundary conditions, a vac- 
uum region of about 10 A was used to avoid interactions 
between a NW and its image. We have considered ax- 
ial lengths (along the NW growth direction) of about 
15.2 and 13.2 A for the [111] and [100] SiC NWs, respec- 
tively. The geometries were optimized using the conju- 
gated gradient scheme, within a force convergence crite- 
rion of 0.05 eV/A. The Brillouin zone was sampled by 
using 2 special k points for 3C-SiC bulk, and 1 special 
k point for [100] and [111] SiC NWs. We verified the 
convergence of our total-energy results with respect to 
the number of special k points using up to 9 k points for 
3C-SiC bulk and 4 k points for [100] and [111] SiC NWs. 

III. RESULTS AND DISCUSSION 

We examined C clusters in 3C-SiC bulk and in 
hydrogen-passivated SiC NWs grown along the [111] and 
[100] directions We considered SiC NWs with di- 
ameter of about 10 A, constructed from the 3C-SiC 
structure£i22 For [100] SiC NWs, two kinds of wires were 
studied: carbon-coated [100] SiC NW and silicon-coated 
[100] SiC NW, whose surfaces are, respectively, carbon 
and silicon terminated. Figure [T] presents the cross- 
section view of the structural models of carbon-coated 
[100] SiC NW [Fig. [Ha)], silicon-coated [100] SiC NW 
[Fig. []Tb)] and [111] SiC NW [Fig. []Tc)]. Note that the 
[111] SiC NW surface has the same number of Si and C 
atoms. 

Eight C defects were investigated 3 - 3 - in the present work 
(see Fig. 02): C Si (carbon antisite), C - C ((100) split of C 
interstitial + C site), C - Csi ((100) split of C interstitial 
+ C antisite), C - Si ((100) split of C interstitial + Si 
site), (C - C)i, 2(C - C) c , [(C - C) Si + (C - C) c ], and 
2(C - C) Si . 

In order to examine the energetic stability of the C 
defects in the 3C-SiC bulk and in the [100] and [111] SiC 
NWs, we used the grand canonical potential at T = K 
(as suggested in the Refs. [H and 1351) 

n = E tot -J2NiPi- (1) 

i 

Here E to t is the total energy of the considered structure, 
fii is the chemical potential of atomic specie i (i — Si, C 
or H) and Ni is the number of atoms i in the system. In 
the thermodynamic equilibrium, 36 

/i Si + ^c=^ k (2) 

and 

/4 u c k = /4 ulk + Mc' lk - AfT(SiC), (3) 




FIG. 1. Cross-section view of the structural models of (a) C- 
coated [100] SiC NW, (b) Si-coated [100] SiC NW and (c) [111] 
SiC NW. Carbon, silicon and hydrogen atoms are represented 
by filled, empty and small-empty circles, respectively. 

where Ai7(SiC) is the formation heat of the SiC bulk. 
Employing diamond structure 3 - 7 - for bulk phases of Si 
and C, and zincblend structur o 37 ' 38 for SiC, we found 
AiJ(SiC) — 0.73 eV, which is in good agreement with 
the experimental value Aff(SiC) = 0.68 e\l£2- 

The upper limits of /xsi and nc are ^g" lk and /iQ Ulk , 
respectively. Hence, the potential fluctuations A/igj = 
MSi — Msi" 114 ano - A/Uc = /j,c — /J.Q Ulk are restricted to 

-AH(SiC)<A MS i,c<0. (4) 

This range defines the C-rich (or Si-poor, where A/ic = 
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FIG. 2. Optimized geometries of the carbon defects investi- 
gated in 3C-SiC bulk and in hydrogen-passivated [100] and 
[111] SiC NWs. (a) C Si , (b) C - C, (c) C - C Si , (d) C - Si, 
(e) (C - C) ; , (f) 2(C - C)o, (g) [(C - C)a + (C - C) c ], and 
(h) 2(C - C) Si . 



and A/xsi = — A-ff(SiC)) and Si-rich (or C-poor, where 
Afi c = -A-ff(SiC) and A^ Si = 0) limits. 

In TableUwe present the formation energies, relative to 
the pristine structures, of C defects in 3C-SiC bulk, and 
in hydrogen-passivated [100] and [111] SiC NWs. The 
two values presented in each column correspond to Si- 
rich (or C-poor) and C-rich limits. For 3C-SiC bulk at 
stoichiometric condition we found that the formation en- 
ergies of C Si , C - C, C - C Si , (C - C)i, 2(C - C) Si and 
[(C - C) Si + (C - C) c ] are 3.49, 7.38, 6.87, 8.79, 9.51 
and 11.55 eV, respectively. These results are in good 
agreement with those obtained by Gali et al. in Ref. HO- 

In 3C-SiC bulk, and in C-coated [100] and [111] SiC 



NWs, the Csi defect [Fig. [2(a)] is the energetically most 
favorable configuration: we find Af2=2.76, 3.05, and 
2.59 eV, respectively. Further formation energy compar- 
ison indicates that, upon presence of interstitial carbon 
atoms, C - Si [Fig. [2(d)] is energetically more stable than 
(C - C)i [Fig. 12(e)], in 3C-SiC bulk, for any C concen- 
tration. We find formation energy differences (AO), be- 
tween C - Si and (C - C) ; , of 1.17, 0.81 and 0.44 eV, 
at Si-rich, stoichiometric and C-rich conditions, respec- 
tively. That is, the formation of (C - C), structures is 
quite unlikely in 3C-SiC. Nevertheless, in both C-coated 
[100] and [111] SiC NWs, (C - C)j is more favorable than 
C - Si at C-rich limit (by 0.21 eV in C-coated [100] SiC 
NW and 0.11 eV in [111] SiC NW). In this case, in con- 
trast with the 3C SiC bulk phase, we may find C - C 
interstitial dimers embeded in SiC NWs [Fig. [2(e)]. How- 
ever, by reducing the concentration of C atoms, we ver- 
ify that C - Si becomes less stable than (C - C)i for 
A/i c > -0.21 and -0.11 eV, in C-coated [100] and [111] 
SiC NWs, respectively. 

Different from C-coated [100] and [111] SiC NWs, 
where Csi is the energetically most favorable configura- 
tion throughout the allowed range for the C chemical 
potential, in Si-coated [100] SiC NW the most stable de- 
fect depends on the C concentration (see Table [J). At 
Si-rich limit C - Si is the most favorable configuration, 
followed by the C - C and (C - C)i defects. However, un- 
der C-rich conditions 2(C - C)si is the most stable defect, 
followed by the C - C Si and [(C - C) Si + (C - C) c ] ag- 
gregates. Figure [3] summarizes our calculated formation 
energy results. The energetically most favorable defect 
in Si-coated [100] SiC NW is C - Si for A/i C < -0.43 eV, 
C - C Si for -0.43 < A/i C < -0.21 eV and 2(C - C) Si for 
A^c > —0.21 eV. At stoichiometric conditions, C - Csi 
is the most stable configuration, followed by C - Si and 
(C - C)j. We find total energy differences, under stoi- 
chiometric conditions, of 0.12 eV between C - Csi and C 
- Si, and of 0.27 eV between C - C si and (C - C);. 

According to Ref. H(J C - Csi is energetically most fa- 
vorable than C - C in 3C-SiC bulk. This result can be 
explained by the larger relaxation required to put two C 
atoms in a C site than to put two C atoms in a Si site. In 
the present work we observe a similar behavior [Fig. 0J: 
under C-rich conditions C - Csi is the most stable com- 
plex with one carbon interstitial in 3C-SiC bulk and in 
[100] and [111] SiC NWs. However, C - C becomes en- 
ergetically more stable than C - Csi for low carbon con- 
centration: A/x c < -0.62 eV (in 3C-SiC bulk), -0.44 eV 
(in C-coated [100] SiC NW), -0.57 eV (in Si-coated SiC 
NW) and -0.51 eV (in [111] SiC NW). At stoichiometric 
conditions, C - Csi is more stable than C - C by 0.52 eV 
(in 3C-SiC bulk), 0.16 eV (in C-coated [100] SiC NW), 
0.42 eV (in Si-coated SiC NW) and 0.29 eV (in [111] SiC 
NW). Our result for 3C-SiC bulk is in close agreement 
with the result reported in Ref. |20j, where C - Csi is 
more favorable than C - C by 0.5 eV, at stoichiometric 
conditions. 

Among the defects with two carbon interstitials, 2(C 
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TABLE I. Relative formation energies (Afi) of carbon defects in 3C-SiC bulk, and hydrogen-passivated [100] and [111] SiC 
NWs. In each column the two values correspond to the Si-rich (or C-poor, where A/^c = — AH (SiC)) and C-rich (A/ic = 0) 
limits (Si-rich/C-rich). The energies are in eV and we used AH(SiC) = 0.73 eV. 

Defect AO. (eV) 





3C-SiC bulk 


C-coated [100] SiC NW 


Si-coated [100] SiC NW 


[111] SiC NW 


Csi 


4.22/2.76 


4.51/3.05 


5.26/3.80 


4.05/2.59 


c - c 


7.75/7.02 


6.78/6.05 


4.94/4.21 


7.09/6.36 


C-Csi 


7.97/5.78 


7.36/5.17 


5.25/3.06 


7.52/5.33 


C - Si 


8.35/7.62 


7.53/6.80 


4.64/3.91 


7.45/6.72 


(C - C), 


9.52/8.06 


8.05/6.59 


5.15/3.69 


8.07/6.61 


2(C - C) c 


12.15/10.69 


10.67/9.21 


6.15/4.69 


10.56/9.10 


[(C - C) Si + (C - C)c] 


10.97/8.05 


9.87/6.95 


6.48/3.56 


9.95/7.03 


2(C - C) si 


13.74/9.36 


12.03/7.65 


6.80/2.42 


12.12/7.74 
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FIG. 3. Relative formation energies (AO,) of the most stable 
C defects in Si-coated [100] SiC NW. The vertical dashed line 
corresponds to the stoichiometric condition. 




0.5 



0.0 



Np-c Si 


c-c - 


: (b) 






c-c - 




\p-c Si J 



Si-rich 



C-rich 



Si-rich 



C-rich 



FIG. 4. Relative formation energies (Afi) of C - C and C - 
Csi defects in (a) 3C-S1C bulk, (b) C-coated [100] SiC NW, 
(c) Si-coated SiC NW and (d) [111] SiC NW. The vertical 
dashed lines indicate the stoichiometric condition. 



- C)si is not expected to occur in 3C-SiC bulk, C-coated 
[100] SiC NW and [111] SiC NW, but is the most stable 
defect, at C-rich limit, in Si-coated [100] SiC NW. Besides 
that, at stoichiometric condition, the formation energy of 
2(C - C) Si is 11.55, 9.84 and 9.93 eV in 3C-SiC bulk, C- 
coated [100] SiC NW and [111] SiC NW, respectively, 
and 4.61 eV in Si-coated [100] SiC NW. One possible 
reason for that is the larger relaxation of the Si-coated 
[100] SiC NW due to its Si-terminated surface. Hence, 
the 2(C - C)s; complex is created with small stress in the 
lattice. This surface effect is expected to disappear for 
larger nanowires. 

We next examine the energy gain upon the formation 
of C clusters through the combination of interstitial C 
atoms. In this case, the energy gain was determined by 
comparing the total energies of the separated systems 
(Ei) and the total of a given (stoichiometrically equiva- 
lent) C cluster (E c ). We define the binding energy of the 
C cluster as E b = E c — Ei . Our calculated binding ener- 
gies are summarized in Table HH The results presented in 
each column correspond to the stoichiometric condition. 
In 3C-SiC bulk or SiC NW with a carbon antisite (Csi 
defect), a carbon interstitial can be captured by C site, 
Si site or C antisite. The last process is more favorable 
than the first (second) one by 4.00 eV (4.60 eV) in 3C- 
SiC bulk, 3.93 eV (4.68 eV) in C-coated [100] SiC NW, 
4.95 eV (4.65 eV) in Si-coated [100] SiC NW, and 3.62 eV 
(3.98 eV) in [111] SiC NW. Considering two carbon inter- 
stitials, the formation of carbon clusters is energetically 
favored in 3C-SiC bulk and in both [100] and [111] SiC 
NW. In fact, the aggregates (C - C)i and 2(C - C) c 
are more stable than two isolated C-C defects by 5.51 
and 2.88 eV in C-coated [100] SiC NW, and by 6.11 and 
3.61 eV in [111] SiC NW. Still, in Si-coated [100] SiC 
NW the energy gain is 13.6 eV to form 2(C - C)si from 
isolated (C - C) and Csi defects, and 3.7 eV to form 2(C 

- C)si from two isolated C - Csi- In 3C-SiC bulk these 
energy gains are 10.2 and 2.2 eV, respectively. Here we 
can infer that, in general, the formation of C clusters is 
quite likely in SiC NWs, in particular for thin Si-coated 
[100] SiC NWs. It is worth to note that, in order to get 



5 



TABLE II. Binding energies (in eV) of C defects in 3C-SiC bulk, and hydrogen-passivated [100] and [111] SiC NWs. In each 
column the values correspond to the stoichiometric condition. 



Cluster 
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Binding 
o-coaieu iiuui 
OlKj IN vv 


energy (eV) 

oi-coaieu iiuui 

d;p mw 
olv^ IN vv 


run <?ip nw 

[111J IN VV 


Cqi+rc - ci rc - 


Cqil 


-4.00 


-3.93 


-4.95 


-3.62 


C Si +[C - Si] ->• [C - 


Csi] 


-4.60 


-4.68 


-4.65 


-3.98 


[c - c ] + [c - q - 


* [(C - C)i] 


-5.98 


-5.51 


-4.73 


-6.11 


[c - c ] + [c - q - 


+ [2(C - C) c ] 


-3.34 


-2.88 


-3.72 


-3.61 


[c - c Si ] + [c - q 


[(C - C) si + (C - C)c] 


-4.75 


-4.27 


-3.71 


-4.66 


[C - C Si ] + [C - C Si 


[2(C - C) Si ] 


-2.20 


-2.69 


-3.70 


-2.92 


2[(C - C)] + 2C Si - 


• [2(C - C) si ] 


-10.20 


-10.55 


-13.60 


-10.16 



a complete picture of the formation of C clusters, the C 
diffusion mechanism is an important issue, however, it is 
beyond the scope of the present work. 
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FIG. 5. (Color online) Electronic band structure and the 
projected density of states (PDOS) of Si-coated [100] SiC NW. 
(a) Pristine NW, and with the presence of defects, (b) C - Si, 
(c) 2(C - C) S i, and (d) C - C Si . In the PDOS diagrams, 
the shaded regions correspond to the total density of states; 
in the electronic band structure diagrams, solid lines (dotted 
lines) represent the spin-up (spin-down) electronic channels. 
The red dashed lines correspond to the position of the Fermi 
level and the zero energy corresponds to the valence band 
maximum. 

Focusing on the electronic properties, we find that sim- 
ilar to the 3C-SiC bulk phase, 20 Csi is an electrically in- 
active defect in SiC NWs. Csi does not introduce states 



within the fundamental band gap. On the other hand, 
the formation of C clusters gives rise to electronic states 
lying within the energy bandgap. Figure [5ja) presents 
the electronic band structure and the projected density 
of states (PDOS) of the pristine Si-coated [100] SiC NW. 
The pristine system exhibits an energy bandgap (at the 
T point) of 1.98 eV, calculated within our DFT-LSDA 
approximation, where the highest occupied (lowest un- 
occupied) states are mostly composed by C 2p (Si 3p) 
orbitals. The formation of C clusters give rise to local- 
ized electronic states within the bandgap, as depicted in 
Figs. EJb)-[5jd) . C - Si defect gives rise to an occupied 
(empty) state at 1.2 eV (1.5 eV) above the valence band 
maximum [Fig [5jbl)]. The dispersionless character of 
those states, along the TX direction (i. e. parallel to the 
NW growth direction) , indicate that those defects states 
are localized around C - Si. Indeed, the PDOS diagram 
[Fig. [5jb2)] indicate that there is a significant contribu- 
tion from the interstitial C atom to the occupied defect 
state, while the nearest neighbor Si atom contributes to 
the formation of the lowest unoccupied defect state. At 
the C-rich limit, the 2(C - C)si defect is the most fa- 
vorable one. Its electronic band structure is depicted in 
Fig. Efcl). We find an occupied state at 1 eV above the 
VBM, and three unoccupied states lying within an en- 
ergy interval of 2.1 and 2.6 eV above the VBM. Those 
defect states are mostly localized around the interstitial 
C atoms [Fig. (5Jc2)]. Finally, for C - Csi (energetically 
most likely at the stoichiometric condition) we find the 
formation of spin-unpaired states within the nearby the 
Fermi level, Figs.[5Jdl) and[5ld2). Those defect states a 
localized around C interstitial atoms. In summary, dif- 
ferent from Csi defects, we verify that the formation of 
C clusters in Si coated [100] SiC NWs gives rise to deep 
states within the bandgap, which may acts as a trap to 
the electronic carriers. 



IV. CONCLUSIONS 

We performed an ab initio investigation of small car- 
bon clusters in SiC bulk and NWs. We examined clusters 
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with carbon interstitials and antisites in 3C-SiC bulk and 
in hydrogen-passivated [100] and [111] SiC NWs. We ob- 
served that the composition of the SiC nanowire surface 
strongly influences the formation of the carbon clusters. 
In fact, Csi is the energetically most favorable configura- 
tion in 3C-SiC bulk, and in C-coated [100] and [111] SiC 
NWs, but is not expected to occur in Si-coated [100] SiC 
NW. The energetically most stable defect in Si-coated 
[100] SiC NW is C - Si at Si-rich limit, C - C Si at sto- 
ichiometric conditions and 2(C - C)si under C-rich con- 
ditions. 

Comparing the total energies of the C - C and C - Csi 
defects (in 3C-S1C bulk and in [100] and [111] SiC NWs), 
we find that C - Csi is more stable than C - C at stoichio- 
metric and C-rich limits. This finding is in accordance 
with the larger relaxation required to put two C atoms in 
a C site than to put two C atoms in a Si site. However, 



we observed that C - C becomes energetically more sta- 
ble than C - Csi for low C concentration. Further total 
energy calculation indicate that the formation of carbon 
clusters in 3C-SiC bulk and in both [100] and [111] SiC 
NWs is energetically favored. Finally, our electronic band 
structure calculations indicate that (i) similar to the 3C 
SiC bulk phase, Csi defects are electrically inactive, while 
(ii) in Si coated [100] NWs the C clusters gives rive to 
deep (localized) levels within the NW bandgap. 
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